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We report on a study of Gilbert damping due to particle-hole pair excitations in conducting 
ferromagnets. We focus on a toy two-band model and on a four-band spherical model which provides 
an approximate description of ferromagnetic (Ga,Mn)As. These models are sufficiently simple that 
disorder-ladder-sum vertex corrections to the long-wavelength spin-spin response function can be 
summed to all orders. An important objective of this study is to assess the reliability of practical 
approximate expressions which can be combined with electronic structure calculations to estimate 
Gilbert damping in more complex systems. 
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I. INTRODUCTION 

The key role of the Gilbert parameter ac in current- 
driveni and precessional^ magnetization reversal has led 
to a renewed interest in this important magnetic ma- 
terial parameter. The theoretical foundations which re- 
late Gilbert damping to the transverse spin-spin response 
function of the ferromagnet have been in place for some 
time&i. It has nevertheless been difficult to predict 
trends as a function of temperature and across mate- 
rials systems, partly because damping depends on the 
strength and nature of the disorder in a manner that re- 
quires a more detailed characterization than is normally 
available. Two groups have recently^ reported success- 
ful applications to transition metal ferromagets of the 
torque- correlation formula^^ for a G - This formula has 
the important advantage that its application requires 
knowledge only of the band structure, including its spin- 
orbit coupling, and of Bloch state lifetimes. The torque- 
correlation formula is physically transparent and can be 
applied with relative case in combination with modern 
spin-dcnsity-functional-thcory-£ (SDFT) electronic struc- 
ture calculations. In this paper we compare the pre- 
dictions of the torque correlation formula with Kubo- 
formula self-consistent-Born-approximation results for 
two different relatively simple model systems, an ar- 
tificial two-band model of a ferromagnet with Rashba 
spin-orbit interactions and a four-band model which cap- 
tures the essential physics of (III,Mn)V ferromagnetic 
semiconductors^. The self-consistent Born approxima- 
tion theory for a G requires that ladder-diagram vertex 
corrections be included in the transverse spin-spin re- 
sponse function. Since the Born approximation is ex- 
act for weak scattering, we can use this comparison to 
assess the reliability of the simpler and more practical 
torque-correlation formula. We conclude that the torque- 
correlation formula is accurate when the Gilbert damping 
is dominated by intra-band excitations of the transition 
metal Fermi sea, but that it can be inaccurate when it is 
dominated by inter-band excitations. 

Our paper is organized as follows. In Section II we ex- 
plain how we evaluate the transverse spin-spin response 



function for simple model ferromagnets. Section III dis- 
cusses our result for the two-band Rashba model while 
Section IV summarizes our findings for the four-band 
(III,Mn)V model. We conclude in Section V with a sum- 
mary of our results and recommended best practices for 
the use of the torque-correlation formula. 



II. GILBERT DAMPING AND TRANSVERSE 
SPIN RESPONSE FUNCTION 

A. Realistic SDFT vs. s-d and p-d models 

We view the two-band s — d and four band p — d mod- 
els studied in this paper as toy models which capture the 
essential features of metallic magnetism in systems that 
are, at least in principle^, more realistically described 
using SDFT. The s — d and p — d models correspond 
to the limit of ab initio SDFT in which i) the majority 
spin d-bands are completely full and the minority spin 
d-bands completely empty, ii) hybridization between s 
or p and e?-bands is relatively weak, and iii) there is ex- 
change coupling between d and s or p moments. In a 
recent paper we have proposed the following expression 
for the Gilbert-damping contribution from particle-hole 
excitations in SDFT bands: 
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where Xxx * s a response-function which describes how 
the quasiparticle bands change in response to a spatially 
smooth variation in magnetization orientation and So is 
the total spin. Specifically, 
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where a and (3 label the x and y transverse spin direc- 
tions and the easy direction for the magnetization is as- 
summed to be the z direction. In Eq.([2]) \i), fi and Uy 
are Kohn-Sham eigenspinors, Fermi factors, and eigenen- 
ergy differences respectively, s a is a spin operator, and 
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Ao(r) is the difference between the majority spin and mi- 
nority spin exchange-correlation potential. In the s — d 
and p — d models Ao(r) is replaced by a phenomeno- 
logical constant, which we denote by Ao below. With 
Ao(r) replaced by a constant jQx reduces to a standard 
spin-response function for non-interacting quasiparticles 
in a possibly spin-dependent random static external po- 
tential. The evaluation of this quantity, and in particu- 
lar the low-frequency limit in which we are interested, is 
non-trivial only because disorder plays an essential role. 



ara formalism, 



X^ P M = -V^j- £ P(i"n,iu n + iw) (3) 



B. Disorder Perturbation Theory 

We start by writing the transverse spin response func- where the minus sign originates from fermionic statistics, 
tion of a disordered metallic ferromagnct in the Matsub- V is the volume of the system and 
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In Eq. (|4|) |ak) is a band eigenstate at momentum k, D is the dimensionality of the system, g(k) = (ak|s x |/3k) 
is the spin-flip matrix element, A Qi/ g(k) is its vertex-corrected counterpart (see below), and 



G a (iu n ,k) = 



Ep — Ek a + i 



sign(w„) 
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We have included disorder within the Born approximation by incorporating a finite lifetime r for the quasiparticles 
and by allowing for vertex corrections at one of the spin vertices. 
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FIG. 1: Dyson equation for the renormalized vertex of the transverse spin-spin response function. The dotted line denotes 
impurity scattering. 



The vertex function in Eq.([4|) obeys the Dyson equation (Fig. (HJ): 

A 0jj g(iw n , iuj n + iu>; k) = s* jj8 (k) + 
f d D k' 

+ J -^u a (k-k')s% ta ,{k,k')G al {iu; n ,k')A a , t p (6) 
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where u a (q) = n a V^(q)(a = 0,x,y,z), n a is the den- 
sity of scatterers, V a (q) is the scattering potential (di- 
mensions: (energy) x (volume)) and the overline stands 
for disorder averagin g 10 ' 11 . Ward's identity requires that 



u a (q) and Tk. Q be related via the Fermi's golden rule: 
— = 2n[ u^k-k'^s^s^JiE^-Ek^,), (7) 

T~aU Jk> , 

a 

where J k = J d D k/{2n) D . In this paper we restrict 
ourselves to spin-independent (a = 0) disorder and 
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spin-dcpcndcnt disorder oriented along the equilibrium- conventional^ integration around the branch cuts of P, 
exchange-field direction(a = zj^. Performing the we obtain 
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t /(e) [P(e + iS,e + iu>) - P(e -iS,e + iw) + P(e -iu),e + iS) - P(e -iv,e- i6)} (8) 



where /(e) is the Fermi function. Next, we perform an 
analytical continuation icu — > u> + it] and take the imag- 
inary part of the resulting retarded response function. 
Assuming low temperatures, this yields 

A 2 

a G = —-^-{Re[P(-id,id)] -Rc \P(id,+iS)]} 

ZTTSq 

= ^L Re (P^ - P R > R ) (9) 
where s = Sq/V, 

pR(A) , R= f A f(/)^(k)G«(0,k)^ iQ (k)Gf^(0,k) 

(10) 

and G R ( A ^ (0, k) is the retarded (advanced) Green's func- 
tion at the Fermi energy. The principal difficulty of 
Eq. ( [9|) resides in solving the Dyson equation for the ver- 
tex function. We first discuss our method of solution in 
general terms before turning in Sections 111 and IV to its 
application to the s — d and p — d models. 

C. Evaluation of Impurity Vertex Corrections for 
Multi-Band Models 

Eq. ( [5]) encodes disorder-induced diffusive correlations 
between itinerant carriers, and is an integral equation 
of considerable complexity. Fortunately, it is possible to 
transform it into a relatively simple algebraic equation, 
provided that the impurity potentials are short-ranged in 
real space. 



Referring back at Eq. ( [5]) it is clear that the solution of 
the Dyson equation would be trivial if the vertex function 
was independent of momentum. That is certainly not the 
case in general, because the matrix elements of the spin 
operators may be momentum dependent. Yet, for short- 
range scattercrs the entire momentum dependence of the 
vertex matrix elements comes from the eigenstates alone: 

< Q ,(k,k') = J2 (ak|m)(m'Kk')C,m< (11) 

m,m' 

This property motivates our solution strategy which 
characterizes the momentum dependence of the vertex 
function by expanding it in terms of the eigenstates of s z 
[s x or s y bases would work equally well): 

A a>/9 (k) = (ak|A|/3k) 

= ]T (ak\m)A mtm ,(m'\pk) (12) 

m,m' 

where \m) is an eigenstate of s z , with eigenvalue m. Plug- 
ging Eqs.f [ITj) and ( [T2|) into Eq.([6|) demonstrates that, 
as expected, A m m ' is independent of momentum. After 
cancelling common factors from both sides of the result- 
ing expression and using <9 q it a (q) = (a = 0, z) we arrive 
at 

A R(A),R_ x ^T TT R(A),R >R(A).R 

1V m,m' — s m,m' m, m':Ll> Iv l,l> K 10 ) 

1,1' 

where 



U%££v = K + ^mmO^(m|ak)G5 A) (0,k)( a k|/}(r|/3k)G«(0,k)(/3k|m') (14) 

I 

Eqs. ( [i"2j) .([i"3" |l and ([HI) provide a solution for the vertex III. GILBERT DAMPING FOR A MAGNETIC 

function that is significantly easier to analyse than the 2DEG 
original Dyson equation. 

The first model we consider is a two-dimensional elec- 
tron gas (2DEG) model with ferromagnetism and Rashba 
spin-orbit interactions. We refer to this as the magnetic 
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2DEG (M2DEG) model. This toy model is almost never 
even approximately realistic*^, but a theoretical study 
of its properties will prove useful in a number of ways. 
First, it is conducive to a fully analytical evaluation of 
the Gilbert damping, which will allow us to precisely un- 
derstand the role of different actors. Second, it enables 
us to explain in simple terms why higher order vertex 
corrections are significant when there is spin-orbit inter- 
action in the band structure. Third, the Gilbert damping 
of a M2DEG has qualitative features similar to those of 
(Ga,Mn)As. 

The band Hamiltonian of the M2DEG model is 
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(15) 



where bk = (— Xk y , Xk Xl Ao), Ao is the difference be- 
tween majority and minority spin exchange-correlation 
potentials, A is the strength of the Rashba SO coupling 



and a = 2s is a vector of Pauli matrices. The correspond- 
ing eigenvalues and eigenstates are 



k 2 

E± M = ~- ± 
2m 
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(16) 
(17) 



-tan 1 {k x /k y ) and 8 



where <J) = 
cos - 1 ( Ao / \J Ag + A 2 k 2 ) arc the spinor angles and 
a = ± is the band index. It follows that 
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a) is a Wigncr function for 
J=l/2 angular momentum^. With these simple spinors, 
the azimuthal integral in Eq. ( 1 14[) can be performed an- 
alytically to obtain 
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where the Kronecker delta reflects the conservation of 
the angular momentum along z, owing to the azimuthal 
symmetry of the problem. In Eq.f I19[) 



cos(0/2) -sin(0/2) 
sin(0/2) cos(0/2) 
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and the retarded and advanced Green's functions are 
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where 
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b k = ^A 2 + A 2 fc 2 , and 7± is (half) 



the golden-rule scattering rate of the band quasiparticlcs. 
In addition, Eq. ( [13)) is readily inverted to yield 
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(22) 



In order to make further progress analytically we as- 
sume that (Ao, Xkp , 7) << Ep = k F /2m. It then follows 
that 7+ ~ 7- = 7 and that 7 = ttN2du° + kN2d\ = 
7o + Jz- Eqs. ([19]) and ([20]) combine to give 



R,R 

-,+:-,+ 
A.R 

-,+:-,+ 
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+,-:+,- 



R,R 
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(70 ~7z) 



— ; cos , 

-b + V 2 



6 + 27 



(23) 



r 



where b ~ ^A 2 fc|, + Ag and cos6> ~ A /b. The first 
and second terms in square brackets in Eq. ( |2"3"]) emerge 
from inter-band transitions (a ^ [3 in Eq. ([19])), while 
the last term stems from intra-band transitions (a = /?). 



Amusingly, U vanishes when the spin-dependent scatter- 
ing rate equals the Coulomb scattering rate ( 7z =70); in 
this particular instance vertex corrections are completely 
absent. On the other hand, when 7z = and b << 7 
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we have U^'+. . — U^'^.^ ~ 1, implying that vertex 

corrections strongly enhance Gilbert damping (recall Eq. 
(I22p), We will discuss the role of vertex corrections more 
fully below. 

A=0.3e F ; A,k F = 0;u°=3u 



T A,R 
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. no vertex corrections 
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FIG. 2: M2DEG: Gilbert damping in the absence of spin- 
orbit coupling. When the intrinsic spin-orbit interaction is 
small, the 1st vertex correction is sufficient for the evalua- 
tion of Gilbert damping, provided that the ferromagnet's ex- 
change splitting is large compared to the lifetime-broadening 
of the quasiparticle energies. For more disordered ferromag- 
nets (Efto < 5 in this figure) higher order vertex corrections 
begin to matter. In either case vertex corrections are signifi- 
cant. In this figure 1/to stands for the scattering rate off spin- 
independent impurities, defined as a two-band average at the 
Fermi energy, and the spin-dependent and spin-independent 
impurity strengths are chosen to satisfy u° = 3u z . 



After evaluating A(k) from Eqs. ([T 
the last step is to compute 



and ( 



p*W,R = f A J^^(k)^ ia (k)G^)(k)Gf(k). (24) 



Since we are assuming that the Fermi energy is the largest 
energy scale, the integrand in Eq. ( [24]) is sharply peaked 
at the Fermi surface, leading to P R ' R ~ 0. In the case of 
spin-independent scatterers {p/ z = — > 7 = 70), tedious 
but straightforward algebra takes us to 



a G (u z = 0) 



A^A 2 (A 2 fcgO(& 2 + Ag + 27 2 ) 
4 So7 o (fe 2 + A 2 ) 2 +4A 2 7 2 ' 



(25) 



Eq. (29) agrees with results published in the recent 
literature^. We note that ac(u z = 0) vanishes in the 
absence of SO interactions, as expected. It is illustrative 
to expand Eq. ( [25f in the b >> 70 regime: 
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which displays intra-band (~ 
70) contributions separately. 



■A 2 ) 70 (6 2 + A 2 )3 

(26) 

7 1 ) and inter-band (~ 
The intra-band damp- 



no vertex corrections 
1st vertex correction 
all vertex corrections 




0.10 
1/(e F x„) 



FIG. 3: M2DEG: Gilbert damping for strong SO interactions 
(XkF — 1.2Ef — 4Ao). In this case higher order vertex cor- 
rections matter (up to 20 %) even at low disorder. This sug- 
gests that higher order vertex corrections will be important 
in real ferromagnetic semiconductors because their intrinsic 
SO interactions are generally stronger than their exchange 
splittings. 
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FIG. 4: M2DEG: Gilbert damping for moderate SO inter- 
actions (Xkp — 0.2Ao). In this case there is a crossover be- 
tween the intra-band dominated and the inter-band domi- 
nated regimes, which gives rise to a non-monotonic depen- 
dence of Gilbert damping on disorder strength. The stronger 
the intrinsic SO relative to the exchange field, the higher the 
value of disorder at which the crossover occurs. This is why 
the damping is monotonically increasing with disorder in Fig. 
(|2| and monotonically decreasing in Fig. (|3]). 



magnetization orientation, the breathing Fermi surface 
effect which produces more damping when the band- 
quasiparticles scatter infrequently because the popula- 
tion distribution moves further from equilibrium. The 
intra-band contribution to damping therefore tends to 
scale with the conductivity. For stronger disorder, 
the inter-band term in which scattering relaxes spin- 
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orientations takes over and ag is proportional to the 
resistivity. Insofar as phonon-scattering can be treated 
as elastic, the Gilbert damping will often show a non- 
monotonic temperature dependence with the intra-band 
mechanism dominating at low-temperatures when the 
conductivity is large and the inter-band mechanism dom- 
inating at high-temperatures when the resistivity is large. 

For completeness, we also present analytic results for 
the case 7 = j z in the b » j z regime: 



Taking 7 = j z for simplicity, we indeed get 



lim olq ~ + 0(7 ) 
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(29) 



a G (u° = 0) 
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where A(l) contains the first vertex correction only, and 
A(oo) includes all vertex corrections. However, the state 
of affairs changes after turning on the intrinsic SO inter- 
action, whereupon Eq. ( I29[) transforms into 



7,66 2 



2A 2 



(3b 2 



A 2 )3 
(27) 

This expression illustrates that spin-orbit (SO) interac- 
tions in the band structure are a necessary condition for 
the intra-band transition contribution to a G - The in- 
tcrband contribution survives in absence of SO as long 
as the disorder potential is spin-dependent. Interband 
scattering is possible for spin-dependent disorder because 
majority and minority spin states on the Fermi surface 
are not orthogonal when their potentials are not identi- 
cal. Note incidentally the contrast between Eq. ( [26]) and 
Eq. (l27|): in the former the inter-band coefficient is most 
suppressed at weak intrinsic SO interaction while in the 
latter it is the intra-band coefficient which gets weakest 
for small A fop. 

More general cases relaxing the (Aq, Xkp,j) << Ep 
assumption must be studied numerically; the results are 
collected in Figs. (HJ, and (gj). Fig (H} highlights 
the inadequacy of completely neglecting vertex correc- 
tions in the limit of weak spin-orbit interaction; the in- 
clusion of the the leading order vertex correction largely 
solves the problem. However, Fig. and ([3]) together 
indicate that higher order vertex corrections are notice- 
able when disorder or spin-orbit coupling is strong. In 
the light of the preceding discussion the monotonic de- 
cay in Fig. ([3]) may appear surprising because the inter- 
band contribution presumably increases with 7. Yet, 
this argument is strictly correct only for weakly spin- 
orbit coupled systems, where the crossover betwen inter- 
band and intra-band dominated regimes occurs at low 
disorder. For strongly spin-orbit coupled systems the 
crossover may take place at a scattering rate that is (i) 
beyond experimental relevance and/or (ii) larger than the 
band-splitting, in which case the inter-band contribution 
behaves much like its intra-band partner, i.e. 0(1/7). 
Non-monotonic behavior is restored when the spin-orbit 
splitting is weaker, as shown in Fig. (2]). 

Finally, our analysis opens an opportunity to quan- 
tify the importance of higher order impurity vertex- 
corrections. Kohno, Shibata and Tatara^ 1 - claim that the 
bare vertex along with the first vertex correction fully 
captures the Gilbert damping of a ferromagnet, provided 
that Aor >> 1. To first order in U the vertex function 
is 



,R(A),R = x tjR(A),R 
m,tn' m,m' ' / 4 m,m':Z,Z ; 
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When Aq << Xkp, both intra-band and inter-band ratios 
show a significant deviation from unityi£, to which they 
converge as A — > 0. In order to understand this behavior, 
let us look back at Eq. (I22p. There, we can formally ex- 
pand the vertex function as A = i X)^Lo wnere the 
n-th order term stems from the n-th vertex correction. 
From Eq. ( [23) we find that when A = 0, U n ~ 0(7") 
and thus n > 2 vertex corrections will not matter for the 
Gilbert damping, which is 0(7)^ when Ep >> 7. In 
contrast, when A / the intra-band term in Eq. ( |23|) 
is no longer zero, and consequently all powers of U con- 
tain 0(7°) and 0(7*) terms. In other words, all vertices 
contribute to 0(1/7) and 0(7) in the Gilbert damping, 
especially if Xkp / Ao is not small. This conclusion should 
prove valid beyond the realm of the M2DEG because it 
relies only on the mantra "intra-band^ 0(1/7); inter- 
band ~ 0(7)" . Our expectation that higher order vertex 
corrections be important in (Ga,Mn)As will be confirmed 
numerically in the next section. 



IV. GILBERT DAMPING FOR (Ga,Mn)As 

(Ga,Mn)As and other (III,Mn)V ferromagnets are like 
transition metals in that their magnetism is carried 
mainly by d-orbitals, but unlike transition metals in that 
neither majority nor minority spin d-orbitals are present 
at the Fermi energy. The orbitals at the Fermi energy 
are very similar to the states near the top of the valence 
band states of the host (III,V) semiconductor, although 
they are of course weakly hybridized with the minority 
and majority spin d-orbitals. For this reason the elec- 
tronic structure of (III,Mn)V ferromagnets is extremely 
simple and can be described reasonably accurately with 
the phenomenological model which we employ in this sec- 
tion. Because the top of the valence band in (III,V) semi- 
conductors is split by spin-orbit interactions, spin-orbit 
coupling plays a dominant role in the bands of these fer- 
romagnets. An important consequence of the strong SO 
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interaction in the band structure is that diffusive vertex 
corrections influence uq significantly at all orders; this 
is the central idea of this section. 

Using a p-d mean-field theory model^ for the ferro- 
magnetic ground state and a four-band spherical models 
for the host semiconductor band structure, Gai_ x Mn x As 
may be described by 



p=0.2nrrf (e F =240meV) ; x=0.08 ; u =3 u' 



H = 



1 

2m 



7i + ^72 ) k 2 - 2 73 (k-s) 2 



+ A s z , (31) 



where s is the spin operator projected onto the 3=3/2 
total angular momentum subspace at the top of the va- 
lence band and {71 = 6.98,72 = 73 = 2.5} are the Lut- 
tingcr parameters for the spherical-band approximation 
to GaAs. In addition, Ao = JpdSNMn is the exchange 
field, J p d = 55mcVnm 3 is the p-d exchange coupling, 



S = 5/2 is the spin of the Mn ions, N 



AIn 



Ax/a J 



is 



the density of Mn ions, and a = 0.565nm is the lattice 
constant of GaAs. 

p=0.6 nm~ 3 (e f =500 meV) ; x=0.04 ; u°=3 u" 
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FIG. 6: GaMnAs: When the spin-orbit splitting is reduced 
(in this case by reducing the hole density to 0.2nm -3 and 
artificially taking 73 = 0.5), the crossover between inter- 
and intra-band dominated regimes produces a non-monotonic 
shape of the Gilbert damping, much like in Fig. (|4}. When 
either 72 or p is made larger or x is reduced, we recover the 
monotonic decay of Fig. ( [5j . 
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FIG. 5: GaMnAs: Higher order vertex corrections make a 
significant contribution to Gilbert damping, due to the promi- 
nent spin-orbit interaction in the band structure of GaAs. 
x is the Mn fraction, and p is the hole concentration that 
determines the Fermi energy Ef- In this figure, the spin- 
independent impurity strength u° was taken to be 3 times 
larger than the magnetic impurity strength u z . 1/tq corre- 
sponds to the scattering rate off Coulomb impurities and is 
evaluated as a four-band average at the Fermi energy. 



The Ao = eigenstates of this model are 

|o,k) =e- ls ^e- lsVe \a) (32) 

where a) is an eigenstate of s z with eigenvalue a. Un- 
fortunately, the analytical form of the Ao / eigen- 
states is unknown. Nevertheless, since the exchange field 
preserves the azimuthal symmetry of the problem, the 
</>-dependence of the full eigenstates |ak) will be iden- 
tical to that of Eq. ( |3"2"|) . This observation leads to 
U m , m ':i,i' oc S m - m i t i-i', which simplifies Eq. ( Q3J. a G 
can be calculated numerically following the steps detailed 



in the previous sections; the results are summarized in 
Figs. ([5]) and ([6]). Note that vertex corrections mod- 
erately increase the damping rate, as in the case of a 
M2DEG model with strong spin-orbit interactions. Fig. 
( O underlines both the importance of higher order ver- 
tex corrections in (Ga,Mn)As and the monotonic decay 
of the damping as a function of scattering rate. The lat- 
ter signals the supremacy of the intra-band contribution 
to damping, accentuated at larger hole concentrations. 
Had the intrinsic spin-orbit interaction been substantially 
weaker—, olq would have traced a non-monotonic curve 
as shown in Fig. ( [5]) . The degree to which the intraband 
breathing Fermi surface model effect dominates depends 
on the details of the band-structure and can be influ- 
enced by corrections to the spherical model which we 
have adopted here to simplify the vertex-correction calcu- 
lation. The close correspondence between Figs. ([5])-([6]) 
and Figs. ( 02K HJ reveals the success of the M2DEG 
as a versatile gateway for realistic models and justifies 
the extensive attention devoted to it in this paper and 
elsewhere. 



V. ASSESSMENT OF THE 
TORQUE-CORRELATION FORMULA 

Thus far we have evaluated the Gilbert damping for 
a M2DEG model and a (Ga,Mn)As model using the 
(bare) spin-flip vertex (a,k\s x \f3,k) and its renormal- 
ized counterpart (a, k| A|/3, k). The vertex corrected re- 
sults are expected to be exact for 1/r small compared 
to the Fermi energy. For practical reasons, state-of-the- 
art band-structure calculations^ forgo impurity vertex 
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corrections altogether and instead employ the torque- 
correlation matrix element, which we shall denote as 
(a, k\K\f3, k) (sec below for an explicit expression). In 
this section we compare damping rates calculated using 
s a p vertices with those calculated using K a ^ vertices. 
We also compare both results with the exact damping 
rates obtained by using A Qi) g. The ensuing discussion 
overlaps with and extends our recent preprint^. 

We shall begin by introducing the following identity^: 

(a,k\s x \(3,k) = i(a, k| [**,«»] |/J,k) 



elements are between essentially degenerate electronic 
states when disorder is treated non-perturbatively^^ 



— {E Ka -E k ^){a,k\s y \P,k) 
-^{a,k\ [H so , S y]\(3,k). 



(33) 



In Eq. ( I33[) we have decomposed the mean-field quasi- 
particlc Hamiltonian into a sum of spin- independent, ex- 
change spin-splitting, and other spin-dependent terms: 
H = Hkin + Hso + H ex , where Hkin is the kinetic (spin- 
independent) part, H ex = Aqs z is the exchange spin- 
splitting term and H so is the piece that contains the in- 
trinsic spin-orbit interaction. The last term on the right 
hand side of Eq. ( I33[) is the torque-correlation matrix 
element used in band structure computations: 




1/(e F x) 



FIG. 7: M2DEG: Comparison of Gilbert damping predicted 
using spin-flip and torque matrix element formulas, as well as 
the exact vertex corrected result. In this figure the intrinsic 
spin-orbit interaction is relatively weak (A&f = 0.05-Ef — 
O.O6A0) and we have taken u z — 0. The torque correla- 
tion formula does not distinguish between spin-dependent and 
spin-independent disorder. 



{a,k\K\/3,k) 



- — (a,k\[H so , S y]\P,k). (34) 
^0 



Eq. ( I33|) allows us to make a few general remarks on 
the relation between the spin-flip and torque-correlation 
matrix elements. For intra-band matrix elements, one 
immediately finds that s% a = K a a and hence the two 
approaches agree. For inter-band matrix elements the 
agreement between ^ and K a p should be nearly iden- 
tical when the first term in the final form of Eq.( |3"3")) 
is small, i.e. when^i (-Ek.a — ^k,/?) << Ao- Since this 
requirement cannot be satisfied in the M2DEG, we ex- 
pect that the inter-band contributions from K and s x 
will always differ significantly in this model. More typi- 
cal models, like the four-band model for (Ga,Mn)As, have 
band crossings at a discrete set of k-points, in the neigh- 
CX 

jj ~ s x a p. The relative weight of 



borhood of which K, 



these crossing points in the overall Gilbert damping de- 
pends on a variety of factors. First, in order to make 
an impact they must be located within a shell of thick- 
ness 1/r around the Fermi surface. Second, the contri- 
bution to damping from those special points must out- 
weigh that from the remaining k-points in the shell; this 
might be the case for instance in materials with weak 
spin-orbit interaction and weak disorder, where the con- 
tribution from the crossing points would go like r (large) 
while the contribution from points far from the cross- 
ings would be ~ 1/r (small). Only if these two con- 
ditions are fulfilled should one expect good agreement 
between the inter-band contribution from spin-flip and 
torque-correlation formulas. When vertex corrections are 
included, of course, the same result should be obtained 
using either form for the matrix clement, since all matrix 



A =0.1e F ; ?ik F =0.5e F ; u =0 




'l/(e F x) 

FIG. 8: M2DEG: Comparison of Gilbert damping predicted 
using spin-flip and torque matrix element formulas, as well as 
the exact vertex corrected result. In this figure the intrinsic 
spin-orbit interaction is relatively strong (A&f = 0.5-Ef = 
5Ao) and we have taken it" = 



In the renaming part of this section we shall focus on a 
more quantitative comparison between the different for- 
mulas. For the M2DEG it is straightforward to evaluate 
etc analytically using K instead of s x and neglecting ver- 
tex corrections; we obtain 



A 



N 2D A Q 
8s 



\ 2 k 2 F A 
b 2 7 



A 2 fc| 
A n 6 



7A0 



T + 



b 2 



(35) 
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where we assumed (7, Afcpi Aq) << c_f. By compar- 
ing Eq. ( [35]) with the exact expression Eq. ( [55]) , we 
find that the intra-band parts are in excellent agreement 
when Ao << Xkp, i.e. when vertex corrections are rela- 
tively unimportant. In contrast, the inter-band parts dif- 
fer markedly regardless of the vertex corrections. These 
trends are captured by Figs. ([7]) and ([H]), which compare 
the Gilbert damping obtained from s x , K and A matrix 
elements. Fig. ( [7J corresponds to the weak spin-orbit 
limit, where it is found that in disordered ferromagnets s x 
may grossly overestimate the Gilbert damping because its 
inter-band contribution does not vanish even as SO tends 
to zero. As explained in Section III, this flaw may be re- 
paired by adding the leading order impurity vertex cor- 
rection. The torque-correlation formula is free from such 
problem because K vanishes identically in absence of SO 
interaction. Thus the main practical advantage of K is 
that it yields a physically sensible result without having 
to resort to vertex corrections. Continuing with Fig. ([7]), 
at weak disorder the intra-band contributions dominate 
and therefore s x and K coincide; even A agrees, because 
for intra-band transitions at weak spin-orbit interaction 
the vertex corrections are unimportant. Fig. ( [5]) cor- 
responds to the strong spin-orbit case. In this case, at 
low disorder s x and K agree well with each other, but 
differ from the exact result because higher order vertex 
corrections alter the intra-band part substantially. For a 
similar reason, neither s x nor K agree with the exact A 
at higher disorder. Based on these model calculations, 
we do not believe that there are any objective grounds to 
prefer either the K torque-correlation or the s x spin-flip 
formula estimate of ac when spin-orbit interactions are 
strong and ag is dominated by inter-band relaxation. A 
precise estimation of etc under these circumstances ap- 
pears to require that the character of disorder, incud- 
ing its spin-dependence, be accounted for reliably and 
that the vertex-correction Dyson equation be accurately 
solved. Carrying out this program remains a challenge 
both because of technical complications in performing 
the calculation for general band structures and because 
disorder may not be sufficiently well characterized. 

Analogous considerations apply for Figs. ( and 
( [TO]) , which show results for the four-band model re- 
lated to (Ga,Mn)As. These figures show results similar 
to those obtained in the strong spin-orbit limit of the 
M2DEG (Fig. [5]). Overall, our study indicates that 
the torque- correlation formula captures the intra-band 
contributions accurately when the vertex corrections are 
unimportant, while it is less reliable for inter-band con- 
tributions unless the predominant inter-band transitions 
connect states that are close in energy. The torque- 
correlation formula has the practical advantage that it 
correctly gives a zero spin relaxation rate when there is 
no spin-orbit coupling in the band structure and spin- 
independent disorder. The damping it captures derives 
entirely from spin-orbit coupling in the bands. It there- 
fore incorrectly predicts, for example, that the damp- 
ing rate vanishes when spin-orbit coupling is absent in 
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FIG. 9: GaMnAs: Comparison of Gilbert damping predicted 
using spin-flip and torque matrix element formulas, as well as 
the exact vertex corrected result, p is the hole concentration 
that determines the Fermi energy Ef and x is the Mn frac- 
tion. Due to the strong intrinsic SO, this figure shows similar 
features as Fig. ([HJ. 



p=0.8nm" 3 (e F =605 meV) ; x=0.04 ; u =0 
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FIG. 10: GaMnAs: Comparison of Gilbert damping pre- 
dicted using spin-flip and torque matrix element formulas, as 
well as the exact vertex corrected result. In relation to Fig. 
( [9]) the effective spin-orbit interaction is stronger, due to a 
larger p and a smaller x. 



the bands and the disorder potential is spin-dependent. 
Nevertheless, assuming that the dominant disorder is 
normally spin-independent, the iGformula may have a 
pragmatic edge over the s^-formula in weakly spin-orbit 
coupled systems. In strongly spin-orbit coupled systems 
there appears to be little advantage of one formula over 
the other. We recommend that inter-band and intra- 
band contributions be evaluated separately when cxq is 
evaluated using the torque-correlation formula. For the 
intra-band contribution the s x and K life-time formulas 
are identical. The model calculations reported here sug- 
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gest that vertex corrections to the intra-band contribu- 
tion do not normally have an overwhelming importance. 
We conclude that ao can be evaluated relatively reliably 
when the intra-band contribution dominates. When the 
inter-band contribution dominates it is important to as- 
sess whether or not the dominant contributions are com- 
ing from bands that are nearby in momentum space, or 
cquivalcntly whether or not the matrix elements which 
contribute originate from pairs of bands that are ener- 
getically spaced by much less than the exchange spin- 
splitting at the same wavevector. If the dominant con- 
tributions are from nearby bands, the damping estimate 
should have the same reliability as the intra-band contri- 
bution. If not, we conclude that the Qg estimate should 
be regarded with caution. 

To summarize, this article describes an evaluation 
of Gilbert damping for two simple models, a two- 
dimensional electron-gas fcrromagnct model with Rashba 
spin-orbit interactions and a four-band model which pro- 
vides an approximate description of (III, Mn)V of fer- 
romagnetic semiconductors. Our results are exact in 
the sense that they combine time-dependent mean field 
theory^ with an impurity ladder-sum to all orders, hence 
giving us leverage to make the following statements. 



First, previously neglected higher order vertex correc- 
tions become quantitatively significant when the intrin- 
sic spin-orbit interaction is larger than the exchange 
splitting. Second, strong intrinsic spin-orbit interaction 
leads to the the supremacy of intra-band contributions in 
(Ga,Mn)As, with the corresponding monotonic decay of 
the Gilbert damping as a function of disorder. Third, the 
spin-torque formalism used in ab-initio calculations of 
the Gilbert damping is quantitatively reliable as long as 
the intra-band contributions dominate and the exchange 
field is weaker than the spin-orbit splitting; if these con- 
ditions are not met, the use of the spin-torque matrix 
element in a life-time approximation formula offers no 
significant improvement over the original spin-flip matrix 
element. 
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